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ABSTRACT 

The heating of the plasma in the solar atmosphere is discussed within both frameworks 
of fluid and kinetic drift wave theory. We show that the basic ingredient necessary for 
the heating is the presence of density gradients in the direction perpendicular to the 
magnetic held vector. Such density gradients are a source of free energy for the excita- 
tion of drift waves. We use only well established basic theory, verified experimentally 
in laboratory plasmas. Two mechanisms of the energy exchange and heating are shown 
to take place simultaneously: one due to the Landau effect in the direction parallel 
to the magnetic held, and another one, stochastic heating, in the perpendicular direc- 
tion. The stochastic heating i) is due to the electrostatic nature of the waves, ii) is 
more effective on ions than on electrons, iii) acts predominantly in the perpendicular 
direction, iv) heats heavy ions more efficiently than lighter ions, and v) may easily 
provide a drift wave heating rate that is orders of magnitude above the value that is 
presently believed to be sufficient for the coronal heating, i.e., ~ 6 ■ 10~ 5 J/(m 3 s) for 
active regions and ~ 8 ■ 10 -6 J/(m 3 s) for coronal holes. This heating acts naturally 
through well known effects that are, however, beyond the current standard models 
and theories. 
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1 INTRODUCTION 

The physical parameters in all three main regions of the 
solar atmosphere (photosphere, chromosphere and corona) 
change both horizontally and with altitude. This by all 
means also includes the collisional frequency, with respect 
to which the solar atmosphere could be, roughly speaking, 
termed as 'strongly collisional' (in the photosphere and in 
the lower part of the chromosphere), 'mildly collisional' (in 
the upper part of the chromosphere), and 'collision-less' (in 
most of the corona and beyond). This classification is not 
absolute as it depends on the ratio of the characteristic 
time t c for certain physical process, and the collisional time 
t u ~ 1/f, where v is the collision frequency. The collisions 
of plasma species include those between plasma particles 
(e-i, e-e, i-i), i.e., Coulomb type collisions (which are typi- 
cal for the corona), and those between plasma particles and 
un-ionized neutrals, e-n and i-n (which are absolutely domi- 
nant in the photosphere and in the lower part of the chromo- 
sphere). In the photosphere, for exampl e, the ion-neutral and 
electron-neutral collision frequencies (|Vranies et al.l 120071 . 
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l2008al : iPandev fc Wardlell2008l ) are about 10 9 and 10 10 Hz, 
respectively, while the ion-ion and electron-ion Coulomb col- 
lision frequencies are around 10 7 and 10 9 Hz, respectively. 
In the corona, the Coulomb collisions occur much less fre- 
quent: v e i ~ 36 Hz and vn ~ 1 Hz (for no = 10 15 m -1 ^ 3 ), 
and u ei ~ 0.4 Hz and m ~ 0.01 Hz (for no = 10 13 m" 1/3 ). 

All plasma modes (with probably one exception only, 
see below) are generally damped by collisions, and therefore 
a strong source is needed whenever one deals with waves in 
the (lower) solar atmosphere. However, such a (local!) source 
sustaining the waves over longer distances or time intervals 
is most often absent, even in models for coronal heating by 
waves. On the other hand, in the practically collision-less 
and very hot corona, Landau damping occurs and damps 
plasma waves very effectively by hot resonant ions, so that 
in this region too, a wave source is required to sustain the 
waves over a longer time interval. 

There is only one mode that is able to survive the drasti- 
cally different (collisional-collisionless) extremes in the dif- 
ferent layers of the solar atmosphere, viz. the drift mode. 
Moreover, this mode is able to benefit (grow) from each of 
these extreme situations. This drift mode has been called the 
'universally growing mode' in the literature. In collisional 
plasma, the drift mode grows due to the electron collisions 
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and this can be described within the two-fluid model. In 
collision-less plasma, however, the mode grows due to the 
electron resonance effect in the presence of a density gradi- 
ent, but this is a purely kinetic effect. The fluid description 
works very well in the lower solar atmosphere, simply due to 
the fact that the plasma is so strongly collisional there (i.e., 
the maxwellization is very effective). On the other hand, 
the kinetic description is preferred in the collision-less en- 
vironment in the corona. Hence, in both extremes we have 
a proper environment and a proper theoretical model for 
studying the strong instabilities of the drift waves. Never- 
theless, the driving mechanism for the waves is the same in 
both cases, namely the presence of a density gradient per- 
pendicular to the ambient magnetic field vector. 

Numerous observations confirm the presence of such 
density irregularities. Yet, the impression is that the pos- 
sible role of the drift wave in the coronal heating problem 
is either overlooked or strongly underestimated (the mode 
is never mentioned in any book dealing with solar plasma). 
This may be partly due to the fact that the drift wave neces- 
sarily implies a multi-component fluid or kinetic description, 
contrary to the widely used single-fluid magnetohydrody- 
namics (MHD) model, within which the drift wave simply 
can not be studied. 

The heating of the corona is one of the long-standing 
puzzles in solar physics and relates to the question of why 
the temperature of the Sun's corona is about 1 to 3 x 10 6 K 
(parts of the corona are even hotter) while the solar surface 
is only around 56 00 K hot. A detailed rev i ew of t he problem 
may b e found in lNarain fc Umschneiderl (ll990T); [ Klimtchukl 
( 2006). According to lNarain fc Umschneiderl (|l990l ) the nec- 
essary heating rates for active regions and coronal holes are 
around 6 • 10~ 5 J/(m 3 s) and 8 • 10 -6 J/(m 3 s), respectively. 
Similar values a nd a detailed analy sis of the problem may 
be found also in I Aschwandenl (|2004l ). T wo theo r ies on coro- 
nal heating, namely wave heating [e.g. ISuzukil {2004)] and 
magnetic reconnection (or nanoflares), have remained as the 
most likely candidates, and many solar physicists believe 
that some combination of these two theories can probably 
explain coronal heating, although 'the details are not yet 
complete'. Both models, however, rely on the continuum ap- 
proximation (MHD) while it is clear that the actual heat- 
ing takes place at length scales much smaller than those 
on which the MHD model is justified. Moreover, it is ev- 
ident that the observed discrepancy between the ion and 
electron temperatures in the corona as well as the observed 
large temperature anisotropy, with a proton perpendicular 
temperature higher than the parallel temperature, are be- 
yond the (sin gle!) fluid model. No te also that according to a 
recent study l|Regnier et al.l [20081 ). the distribution of mag- 
netic null points (only 2% of them are located in the corona 
and 54% in the photosphere) is opposite to what would be 
required for the mechanism that is supposed to heat the 
corona. The heating by waves rather than by reconnection is 
al so supported by the d iagnostic of active regions presented 
in iMilligan et "all (|2005l ). 

A self-consistent heating model must fulfil a lot of re- 
quirements. It must: 1) provide an energy source for the 
extremely high temperature in corona, including 2) a reli- 
able and efficient mechanism for the energy transfer from the 
source to the plasma particles, and 3) this with a required 
heating rate. It should also 4) explain the discrepancy be- 




Figure 1. TRACE photos of numerous loops in active regions, 
with clear density structures. 



tween ion and electron temperatures (typically Ti > T e ), 
5) explain the origin of the large temperature anisotropy 
(T± > Tii) with respect to the direction of the magnetic 
field, particularly for ions, 6) explain the observed larger 
heating of heavier ions, and last but not least, 7) it should 
work everywhere in corona (with well known different heat- 
ing requirements in active and quiet regions). 

In this paper, we present a new coronal heating model 
that can operate in all layers and in all magnetic struc- 
tures of the solar atmosphere and that is able to explain all 
seven requirements given above. Our model represents a new 
paradigm in which a) the energy for driving the drift modes 
and for the heating of corona is already present in corona, 
and, b) this energy is naturally transmitted to the different 
plasma species by well known effects that are, however, be- 
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yond the standardly used models and theories. Moreover, it 
is based on well established, basic theory which has already 
been verified and confirmed by means of laboratory plasma 
experiments. All that is needed for the heating mechanism 
to work is the presence of a density gradient perpendicular 
to the magnetic field. In the solar corona, this may be taken 
rather as a fact than as a hypothesis. 



2 INHOMOGENEOUS SOLAR PLASMA 

Fine density filaments and threads in the solar atmosphere 
have been observed for a long time now, even from ground- 
based observations like those during the eclipse in 1991 
(INovember fc Koutchmy|[l997h . showing a slow radial en- 
largement of plasma structures. On the ot her hand, mea- 
surem ents by Voyager 1 and Voyager 2 show ( Wo o fc Habball 
1 19971 ) that the finest structures in the slow solar wind at 
around 9 Rq are about 3 times finer than those in the 
fast wind. Assuming a radial expansion, these authors con- 
clude that the transverse sizes of these highly elongated 
structures at the Sun are below 1 km. T he contour maps 
presented in iKarowska fc Habball |l99ll) reveal the exis- 
tence of numerous structures of various size. The small- 
est filamentary s tructures of the order of 1 km have been 
discussed also in IWool dl996h . Recent Hinode observations 
(|De Pointieu et al.ll2007 ) confirm that the solar atmosphere 
is a highly structured and inhomogeneous system and re- 
vealed radially spreading grass-like density filaments of var- 
ious size pervading th e whole domain. A very r ecent three- 
dimensional analysis l|Aschwanden et all l2008ah of coronal 
loops reveals short-scale density irregularities within each 
loop separately. Particularly clear imag es of coronal loops 
with de nsity struct ures may be s e en in lAschwanden et aTJ 
l|2008bl l and also in lWarren et alj (|20Qgt ). In Fig. 1 we give 
some TRACE images showing many coronal loops in active 
regions. Such active r egions do not cover the whole corona, 
yet the analysis from lAschwanden! (|200ll ) reveals that they 
require around 82% of the total energy needed to heat the 
entire corona. 

The presence of plasma density irregularities through- 
out the corona implies a plasma that is not in thermody- 
namic equilibrium. In other words, it reveals the presence 
of free energy in the system. In all the examples mentioned 
above, such density irregularities are as rule associated with 
the magnetic field, creating the ideal environment for drift 
waves. The purpose of this work is to show that the energy 
stored in these density gradients may drive drift waves on 
massive scales. We shall point out some basic features of the 
drift wave instability, apply them throughout the solar atmo- 
spheric plasma and compute the appropriate growth rates. 
It will be shown that the growing drift wave and its subse- 
quent interaction with plasma particles may yield the long 
searched solution of the problem of the heating of the solar 
atmosphere. The model implies that the direct energy supply 
for the heating comes from within the corona itself, though 
still maintained and replenished by some mechanisms from 
below the photosphere. Those include a continuous restruc- 
turing of the magnetic field implying the consequent similar 
changes of the plasma density (due to the frozen-in condi- 
tion), and also t he observed inflow o f the plasma along the 
magnetic loops l|Schriiver et al.|[l999h . To some extent, this 



looks similar to the currently accepted scenarios mentioned 
above, where the magnetic field plays an essential role and 
is assumed as given. However, in this new approach the dis- 
sipation of these drift waves is easy to explain in our kinetic 
model that works on the (very small) length scales at which 
the actual dissipation takes place. 

The observable characteristic dimensions of the den- 
sity irregularities are limited by the available resolution of 
the instruments (about 0.5 arcsecond in the example from 
Fig. 1, that is below 400 km on the Sun). Nevertheless, 
even extremely short, meter-size scales can not b e excluded, 
especially in corona l|Vranies fc Poedtsl l2008al ). This can 
be seen by calculatin g the perpendicular ion diffusion co- 
efficient l|Chenl fl988f ) for a coronal environment: Dj_,j ~ 
KTjfj I (mjCllj) oc rn~J 2 , where j = i, e. Taking Bo = 10~ 2 T, 
n = 10 15 m~ 3 , and T e = T % = 10 6 K, for ions we obtain 
Dx,i = 0.01 m 2 /s. The diffusion velocity in the direction 
of the given density gradient is (|Vranies fc Poedtsl 120061 ) 
DxjVno/no [see also Eq. © further in the text]. Taking 
the inhomogeneity scale-length L n = [{dno/dx)/na\~ l — 10, 
10 2 , and 10 5 m, where ^-denotes the direction perpendicular 
to the magnetic field vector, we obtain for the ion diffusion 
velocities, respectively, 10 -3 , 10 -4 , an 10 -7 m/s only. There- 
fore, even very short density inhomogeneities can last long 
enough to support relatively high frequency drift instabili- 
ties. Hence, in dealing with the drift wave, we may operate 
with the density inhomogeneity scale lengths that have any 
value ranging from one meter up to thousands of kilometers 
in the case of coronal plumes. 



3 DRIFT WAVE WITHIN FLUID THEORY 

Assuming a partially /weakly ionized and collisional plasma, 
like in the photosphere and chromosphere, it is justified to 
employ the fluid model. In such an environment, the kinetic 
Landau damping is not expected to play any significant role 
as long as the ion mean free path is below the wavelength. 
The strong-weak La ndau damping t ransition has been ex- 
perimentally verified (|D'Angelolll968l ) to be at the threshold 
uj ~ Vi. The momentum equations that we use for electrons 
and ions can be written as 



9vi — _, 

+ (Vi ■ V)Vi 



- uTiVlli — V • 7Ti 



rrnniViVi, 



and 
= 



+ VeX Bo) - kT e Vn e 



(1) 



(2) 



respectively. Here, Vi 



and v e = v m + v ei . The ion mo- 



mentum change due to collisions with electrons is neglected 
in view of the mass difference. This is justified as long as the 
ion dynamics in the direction of the magnetic field vector is 
neglig ible. Otherwise it may modify the instability thresh- 
old (jVranies fc Poed ts 2008b). The shape of the equations 
reveals that we are dealing with electrostatic perturbations, 
the hot ion effects are included through the pressure and the 
gyro-viscosity stress tensor terms, while the left-hand side of 
the electron momentum equation is omitted, implying per- 
turbations with phase speed and perturbed velocity both 
much below the electron thermal velocity. The dynamics of 
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Here and further we have taken Bo = BqS z , and v^i = 
KTi I mi . 

The ion perpendicular velocity obtained from Eq. JT} is 
described by the following recurrent formula 

1 _ „ , , v%i _ Vj_n, Vi V±(j> 
— e z x V±(j>+ — e z x — — — — 

-DO "i Mi 1 ii -DO 



^if T i V±7lj ^ V 

he z x 

Vi Vj_ • 7Tj 1 Vi ( d 

fi; erii B t 



+ »i-V e z X(?,j 



(3) 



Here, aj = l/(l + i/^/flf). The velocity can be calculated up 
to small terms of any order using the drift approximation 
|d/dt|<fi<. 



Figure 2. A drift wave in cylindric geometry with poloidal wave 
number m = 2. 



the ne utrals may also be includ ed. However, as shown else- 
where l|Vranies fc Po edts 2009), this usually yields small or 
negligible corrections, and such a model works well as long 
as the ion sound response is negligible. Otherwise, there is an 
instability threshold that is modified whe n the dynamics of 
the ne utrals is self-consistently included l|Vranies fc Poedtsl 
2008b). A typical geometry of the drift wave in cylindric 
coordinates is presented in Fig. 2. Here, as an example we 
consider a wave propagating in the poloidal and axial direc- 
tions, with the poloidal mode number m = 2, in a plasma 
with a radially dependent equilibrium density. The actual 
wave fronts are twisted around the axis and have an r- 
dependent amplitude which reaches a maximum in the area 
of the largest density gradient. The mode behavior in the 
presence of both ra dial and axial equi l ibrium density gradi- 
ents is discussed in IVranies fc Poedtsl (|2007f ) . 

In some regions, the drift mode may become electro- 
magnetic, provided the plasma /3 exceeds the electrostatic 
limit. This yields a coupling of the drift and the k inetic 
Alfven wave (|Kadomtsevl 19651 ; I Vranjes fc Poedtsl2006l ) . The 
first experimental observation of such a coupled mod e in a 
hot-ion plasma is discussed in INishida fc IshiH |T974). The 
drift-wave part of the observed coupled modes, is strongly 
growing, and the maximum growth rate u)i ~ 0.1 u) r , where 
uj t is the drift wave frequenc y. In the case of the low er so- 
lar atmosphere, as shown in IVranies fc Poedtsl (|2006l ). the 
Alfven part of the mode is of no interest as it is always 
damped due to the collisions, and it will not be discussed 
here. Note that the usual gas viscosity is as a rule negligi- 
ble even for a relatively high density environment like the 
p hotosphere and chrom osphere. More details are available 
in lVranies etall i|2008al ). 

The ion gyro-visc ous component s of the stress tensor 
that we need here are (|Weilandll2OO0l ) 



riiKTi ( dv 



2Qi 



i).r 



dviy 
dy 



71% KjTi 

~20~ 



OViy 

dx 



dy 



3.1 Ion diamagnetic current effects 

The ion gyro-viscosity is usually overlooked in the litera- 
ture and, in particular, in the domain of solar plasma, and 
this even in studies using the multi-component plasma the- 
ory. For the solar plasma case, this may have very important 
consequences because ion and electron temperatures are typ- 
ically of the same order, so that a cold ion approximation 
and the consequent neglect of the ion gyro- viscosity can not 

be justified. 

As a matter of fact, it is well known (|Vranies fc Poedtsl 
120061 ; IWeiland|[2000l ) (yet standardly disregarded) that part 
of the ion gyro-viscosity contributes to the cancelation of 
some terms in the ion continuity equation. Clearly, this im- 
plies that, if the gyro- viscosity terms are omitted in the 
derivation, the resulting ion equations contain terms that 
can not possibly be there, yielding some false physical ef- 
fects. The cancelation of terms appears when Eq. ([3]) is set 
into the ion continuity dni/dt+X7 ±-(niVi±)+V z -(riiVi z ) = 0, 
and it involves, respectively, the diamagnetic and polariza- 
tion drifts from one side 



e z x v i± , 



and the stress tensor drift term v ni — e z x Vi • 7r \ / '(miUifli) 
from the other side. 

It is seen that, as long as the magnetic field is homoge- 
neous, we have 



[riiVti 



0. 



(4) 



describing a well known fundamental property. This is due to 
the fact that the diamagnetic drift is a fluid effect and not a 
particle drift, and therefore it can not contribute to the flux 
in the continuity equation. It appears due to the gyration of 
ion particles in the presence of density gradient, and without 
any macroscopic motion of the ion guiding center. However, 
if an inappropriate linearization is done, it gives rise to the 
terms (k ■ Vn)ria in the continuity equation, that provides 
a false source of the current-driven instability that may be 
seen in the literature. 

The mentioned cancelation of terms is due to the con- 
vective derivative part 



(Hi ■ V)e z x v i± 



(5) 



in the p olarization drift v V j. The procedure is des cribed in 
detail in IVranies fc Poedtsl (|2006h ; IWeilandl (|200GT ). Within 
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the approximation of small gradients of the equilibrium 
quantities, the last vt± in the convective derivative ([5]) con- 
tains only the leading order perturbed drifts from Eq. Q. 
On the other hand, the first Vi in ([5]) can only be the equi- 
librium ion diamagnetic drift. This is then to be used in the 
term Vj_ • (riiVpi) in the continuity equation. 
The stress tensor drift term yields 



Vi ■ (m?^) = —pi\7±n i0 ■ V±ViA 



■ n i0 p?ViVj 



(6) 



Here, pi — v T i/0,i. Within the second-order small terms 
approximation, the first term on the right-hand side in this 
expression cancels out exactly with the contribution of the 
above discussed convective derivative in the polarization 
drift X7± ■ [n(vi ■ V)e z x Vi±]. The cancelation of terms is ex- 
act and valid for any plasma. Moreover, the results obtained 
from this formal fluid theory can easily be obtained by us - 
ing the kinetic theory as well [cf. IVranies fc Poedtsl |2006); 
IWeilandl(|2000h ]. However, if the derivation is performed in- 
correctly, by simply ignoring the stress tensor contribution, 
then the resulting equations contain extra terms originating 
from the convective derivative in the ion polarization drift 
For perturbations of the form ~ exp(— iuot + ik y y + ik z z 
this implies the presence of terms like k ■ i7«, which in real- 
ity cancel out exactly. A drastic example where a false 'new 
instability' is obta i ned d ue to these terms may be seen in 
iMecheri fc Marschl (120081 ). 



3.2 Dissipative instability 

In the limit of perturbations with a parallel (to the mag- 
netic field vector) phase velocity that is considerably larger 
than the sound speed ui/ k z ^> c 3 = (^T e /mj) 1 ^ 2 , the ion 
continuity equation reads (|Vranies fc Poedtsll200d ): 



d ( n 
dt 



(^)+ — 

V n i0 ) B 



e z x V I 



riio 



2„ 2 ( e<j)i n n \ 2 a 2 /e<^i nn 
-ViPi V_l — - H Pi — Vj_ ( H 



V KTi UiO 



1 dt V nTi riio 



Q,i dt 



Bq fii riio 



0. 



(7) 



Here, the discussion from the previous section is included 
self-consistently. For the assumed shape of perturbations, 
Eq. Q yields: 



nn _ u*j + ojp 2 ky(l + pjky) + ivjp\k 2 y e0i 
n i0 lo[1 + p\k%{\ + p\k%)\ + iv l p\k% uTi' 



(8) 



U T i ILiO ^_ I — * i /i 

cj*i = ky — , \/±nio = — e x n i0 — —e x dnio/dx. 

Hi UiO 

The electron perpendicular and parallel velocities are ob- 
tained from Eq. (J2J: 



1 



i + ui n /m 



Or 



Bo Bo 



T -e z x 

He n e 



V ez l 



ik z v Te / e(j>i n e 



Ven V K± e n e , 

From the electron continuity we then obtain 



(9) 
(10) 



n e i _ uj, e + iD p + iD z efii 
n e Q to + iD p + iD z nT e 

Here D p = vjt 2 y p\, D z = k 2 z v 2 e /v e , p e = v Te /U e , and 



(11) 



kii 



The term D p describes the usually neglected effects of elec- 
tron collisions in the perpendicular direction. In most cases, 
in view of the small inertia, these electron collision effects 
are included from the electron parallel momen tum. This is 
justified provided that (|Vranies fc Poedt3l2006l ) 



k 2 z V^,l{k 2 yV 2 ) > 1. 



(12) 



In solar plasmas, for perturbations with an almost arbitrar- 
ily small parallel wave number this condition may not always 
be satisfied. This means that collisions must be taken into 
account in the perpendicular dynamics. 

In the absence of an equilibrium electric field and for 
quasi-neutral perturbations, the dispersion equation of the 
drift mode in collisional solar plasma is: 

) Te U ti + Up 2 k 2 y(l + P 2 k 2 y) + W t p 2 k 2 y 

T % w[l + p?fc»(l + pikl)\+Wip\k% 

I u*e+i(D P + Dz) _ n m . 
+ lo + i(D p + D Z ) (i6) 

In various limits, Eq. (113[) yields different dispersion equa- 
tions for the drift wave known from the literature. 

The most simple case is for collision-less plasma, when 
the right-hand side in Eq. (|13p reduces to 1. Assuming, in 
addition, that 1 S> p 2 k 2 , we have 

LU r — W» e /(1 + Plk 2 y ), p a =C s /Qi, C*=KT e /mi. (14) 

In the case of cold ions and for D z S> D p , we obtain from 
Eq. (|13p the same oscillation frequency as in the previous 
e xample . But the growth rate due to electron collisions is 



l|Vranies fc Poedtsll2008bl : IWeilandll2000h 

2 2 ; 2 
V e LU r p 8 k y 



(15) 



As mentioned earlier, this describes the drift mode growing 
in the presence of electron collisions. The full Eq. (|13p can be 
solved numerically, taking relevan t parameters for the sola r 
atmosphere at various altitudes (|Vranies fc Poedtsl [2006). 
The result represents an interplay between the ion collisions 
that damp the mode, and the electron collisions that pro- 
duce the growth rate approximately given by Eq. (|15p . This 
may have a great importance in terms of plasma heating be- 
cause the wave is always present in inhomogeneous plasmas, 
and it grows on the account of collisions of light species, 
while in the same time the wave energy is continuously ab- 
sorbed by the collisions of heavy species. The process con- 
tinues till it runs out of energy. In other words, it continues 
as long as the source (i.e., the density gradient) is present. 



4 KINETIC INSTABILITY OF DRIFT WAVE 

In most of the corona the collisions are not expected to 
play an important role and consequently the drift dissipa- 
tive instability may not be of much importance. In addition, 
regarding the problem of coronal heating as a background 
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and motivation for any wave analysis, in the case of so lim- 
ited collisions an efficient mechanism for transfer of energy 
from the wave to plasma is missing. However, this is not so 
within the frame of the drift wave kinetic theory. In that 
case, the process develops as follows: the interaction of the 
wave and electrons is destabilizing and the mode grows due 
to Cherenkov-type interaction (in the presence of the den- 
sity gradient) , while in the same time its energy is absorbed 
by ions d ue to Landau damping . This may be seen from 
IWeilandl (|200Ch : llchimarul | 19801 ) where, under the condi- 
tions 



k z v T i <C u> -C k z v Te , lo 



the wave properties are described by the frequency 

LJ*iAa(bj) 

Wr ~ 1 - Ao(6i) + Ti/T e + k%X 
and the corresponding growth rate 



2 ' 
di 



(17) 



7 ; 



G) 



1/2 



|w,i|Ao(6») 



T e \k z \v Te 



■ exp[-u3 r /(k K v Te )] 



V T i 



■ exp[-u) r /(k z v Ti )] 



(18) 



Here, 



Ao(bi) = Io(bi) exp(— bi), bi = k y pl, X di = v Ti /Lo pi , 

and 7o is the modified Bessel function of the first kind and of 
the order 0. The first condition in 1)160 refers to the use of ap- 
proximative expressions for the plasma dispersion function, 
while the second condition implies a strongly magnetized 
plasma, and the third implies that the acoustic part in the 
dispersion equation can be omitted. Eqs. (|17[) and Q8) , as 
well as Eq. (|13[) . are obtained using a local approximation. 
This implies that the characteristic length for the change of 
the mode amplitude in the direction of the density gradi- 
ent, is much larger than the perpendicular wavelength, i.e., 
(d/dx)^ 1 ~ L n S> X y . Otherwise, an eigen-mode ana l ysis is 
to be used, see for more details in IVranies fc Poedts! (120041 ') 
and IVranies fc Poedts! (|2007T l. Note that in the appropriate 
limits {kl\ 2 M < 1, bi < 1) Eqs. {2} and {T7J) coincide. In 
Eq. (|18[) , the last term in the square bracket is positive. It 
describes the damping on ions and introduces a threshold in 
the mode instability that is of importance only in plasmas 
with hot ions, yet negligible as long as ui r /k z S> v T i (hence 
the reason for interest in small k z ). The necessary condi- 
tion for the instability follows from the first term where we 
must have uj r < to* e , that is as a rule easily satisfied. Equa- 
tions (|14|) and (|17[) reveal the presence of the energy source 
already in the real part of the frequency u> r oc Vino, while 
details of its growth due to the same source are described 
by Eqs. (fT5)l and (fl8)) . 

As an example, in Fig. 3 the growth rate (|18[) is calcu- 
lated for an electron-proton plasma in terms of the parallel 
wavelength \ z by taking B = 10" 2 T, n = 10 15 m~ 3 , 
Ln = [(dno /dx) /no] -1 = s ■ 100 m, and for X y — 0.1, 0.3, 
0.5, and 1 m. Because k y 3> k z , the corresponding wave 
frequencies are practically constant and have approximate 
values 186, 210, 254, and 307 Hz, respectively (this assum- 
ing s = l, but see below). In all these cases, we have very 
strongly growing drift modes. Note that in all four cases 



Figure 3. The growth rate JTHJ normalized to the wave frequency 
u) r in terms of the parallel wavelength, for several perpendicular 
wavelengths X y , and for the density scale-length L n = s ■ 100 m, 
s G (0.1, 10 3 ). 



uj r < w* e , and the growth rate changes sign because of the 
varying (with X z ) contribution of the ion part in Eq. (|18[) . 
Observe the extremely fast growth for short perpendicular 
wavelengths X y . The growth rate increases with A z and may 
easily become much larger than ui r . Yet, strictly speaking, 
in that case the assumption of smallness of the imaginary 
part with respect to the real part, used in the expansion of 
the plasma dispersion function, is violated, and the problem 
must be treated numerically. In the examples presented in 
Fig. 3, the conditions used in obtaining Eq. 1)180 are formally 
well satisfied, i.e., u) T <C Qi and in the given range of k z , we 
have k z v T i < 0J r < k z v Te , and also the sound branch is al- 
ways far below the drift branch. For the given parameters, 
the Debye length Ad is around 1 mm only, and even shorter 
values of X y are permissible. The local approximation is also 
well satisfied and the mode growth is expected to take place 
throughout the density gradient. For the given density, the 
plasma-/? is around 0.6m e /irii. 

A simple way to demonstrate that a similar mode be- 
havior can take place at various inhomogeneity scale-lengths 
L n (in other words at various places in the corona) , is to keep 
the ratio X z /L n fixed. Thus, with the same parameters as 
above we set L„i = s ■ L n , X z i = s • X z , where s takes values, 
e.g., between 0.1 and 10 3 . It can easily be shown that the 
graphs from Fig. 3 remain exactly the same. In other words, 
the ratio 'y/co r remains unchanged, although both 7 and u) r 
are shifted towards lower values. For example, for s = 10 3 , 
i.e., L„i = 10 5 m, X z i — 10 r m, and taking X y = 0.5 m, 
we have 7/oj r = 0.26. This is the same as the value of the 
line labeled c in Fig. 3, for X z = 10 km (and correspondingly 
L n = 100 m). Yet, now 7 = 0.07 Hz, and u r = 0.25 Hz. Such 
a variation of s may be used to describe the natural change 
of the radial density gradient when we move along a mag- 
netic flux tube, that may appear due to the diverging tube 
geometry. So the mode growth occurs everywhere along the 
given flux tube. Hence, as long as the conditions 1)160 used 
in the derivations are satisfied, the results in Fig. 3 are valid 
for any s in the range 0.1 — 10 3 . 
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Figure 4. The growth rate H18H normalized to u) r in terms of the 
perpendicular wavelength, for several values of A 2 . 




Figure 6. The growth rate ( 1181 normalized to the wave frequency 
uj r in terms of the perpendicular wavelength X y and the density 
scale-length L n , for A z = s ■ 2 • 10 4 m, s g (0.1, 10 3 ). 



from Fig. 6, where the growth rate is plotted in terms of 



Figure 5. The growth rate l|18|l normalized to u> r in terms of the 
density scale length L n for several values of A z , s g (0.1, 10 3 ). 



Next, we check the mode behavior with respect to the 



perpendicular wavelength X y . We fix L n 



100 m, keep 



the other parameters the same as above, and calculate the 
growth rate for A z /s = 7, 10, and 15 km. The result is 
presented in Fig. 4. The graphs and the damping in the 
short X y range (for \ z — 7 m) are in agreement with Fig. 3. 

In Fig. 5 we present the mode growth rate in terms of 
L n , for the same plasma parameters as above. The graphs 
remain unchanged for any s € (0.1, 10 3 ). The frequency is 
507 Hz (at L n = 50 m), and 254 Hz (at L n = 170 m), for 
s = 1. In the case s = 10 3 , the frequency is 0.5 Hz (where 
now L n = 50 km), and 0.15 Hz (at L n — 170 km). 

A similar quantitative analysis may be preformed for 
densities several orders of magnitude below the used value, 
and for at least one order of magnitude higher density as 
well, resulting in the same wave behavior, with only some 
shift in the wave frequencies and wavelengths. This is seen 



A y and s ■ L„ by taking no = 10 



Bo 



10-' T, and 



A z = s ■ 2 • 10 m. The smaller values of no and Bo imply the 
possible application of the model to higher altitudes in the 
corona. It is seen that the perpendicular wavelength is just 
shifted up for one order of magnitude, and the growth rate 
decreases with X y and L n just like in Figs. 4 and 5. 

On the other hand, for densities that are two or more 
orders of magnitude larger than the earlier used value 
10 15 m -3 , and thus for a larger plasma-/?, a coupling with 
electromagnetic Alfven-type perturbatio ns may take place 
(lHasegawa I1975T 1. As shown elsewhere l|Vranies fc Poedtsl 
2006), the drift mode behavior will remain similar even in 
that case. The frequency becomes slightly reduced and a 
part of the wave energy is spent on the coupling with this 
additional Alfven mode that appears to be always damped 
by collisions. This coupling is de s cribed by the following d is- 
persion equation 



> coupling is described by the following di 
(|Weilandll2000l ; IVranies fc Poedtsfe oOG) 



Ul — lj 2 (lj„ b + LU ti ) + Lu[bJ* e 0J,i — k z c a /(l + kyPi ) 

- klk 2 z c a (p 2 i + p 2 s )] + k 2 z c 2 a Lu te /(l + k 2 y p 2 ) = 0. (19) 

Here, c? a = Bo / {ponoirii). The nature of coupling is best seen 
for negligible ion thermal effects when we have (lj— Lu„e){oJ 2 — 
k z Ca) — Ldk 2 Cakyp 2 = 0. For sufficiently small k 2 p 2 , the two 
modes propagate practically almost independently. 

In view of all these results, we stress again an important 
difference of the drift wave as compared to other plasma 
modes. The mode frequency given by Eq. (|17[) . already by 
its form implies the presence of an energy source stored in 
the density gradient. This free energy is then responsible 
for its growth, either due to kinetic or fluid effects, and the 
details of the growth are described by Eqs. (|15[1 and (|18[) . 
Compare this with, e.g., an Alfven u> — kc a , or a sound mode 
lj — kc 3 , where the given dispersion equations only allow for 
the possibility for the plasma to support these modes, on 
the condition that an additional energy source is provided. 

As for the experimental verification of such a very 
strongly growing drift-wave instability, one (out of many) 
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wa ve-electric field Eg cos(k z z —u r t), is approximately given 
bv lBittencourtl (jl995T ) 



Figure 7. Velocity of electrons (1201 ) accelerated by the drift-wave 
electric field with <f> = 80 V, for three different starting velocities. 



may be found e.g., in iBrossier et al.l (|l97lh for a similar 
almost collision-less hydrogen plasma with T e = 14 eV, 
Ti = 2 eV, B < 0.35 T, and the density n = 10 16 m -3 . 
Note that in the experiment the maximum observed linear 
growth-rate was very high: u>i ~ ui r — 7 • 10 5 Hz. 



4.1 Electron acceleration by parallel wave-electric 
field 

The electrostatic drift mode presented in Figs. 3-6 im- 
plies a time- varying electric field, whose parallel component 
\X7 z (j)\/E d = \k z (f>\/Ed, normalized to the Dreicer runaway 
electric field, for A z = 10 km and (f> — 60, 80 V/m, is 3.1 and 
4.1, respectively. Here, we use X y = 0.5 m, Bo = 10~ 2 T, 
n = 10 15 m~ 3 , L n = [(dn / dx) / no]' 1 = 100 m, and the 
reason for the given values of <f> will be give n in the forth - 
coming sections. The Dreicer electric field is (|Dreicerlll959r i 
E d = eL e i/(4n£ \j). Here, L ei = log(A d /&o) is the Coulomb 
logarithm, \ d — \ de X di /(\ de + A^) 1 ' 2 is the plasma Debye 
radius, and bo = \e 2 /\12-KEgK (T e -\-Ti)\ is the impact parame- 
ter for electron- ion collisions IVranies et al.l (|2008bh . For the 
given temperature of one million K, we have L e i = 20.1, 
\ d = 0.0015 m, and the Dreicer field is 0.012 V/m. Hence, 
the parallel wave field exceeds the Dreicer field so that the 
bulk plasma species (primarily electrons) can be acceler- 
ated/decelerated by the wave in the parallel direction. The 
acceleration is more effective on particles that are already 
more energetic, resulting in a distribution function consid- 
erably different from a maxwellian. This may be one of the 
reasons behind the observed kappa-distribution in the outer 
solar atmosphere and in the solar wind. 

For s = 10 3 and for the two values of cf> mentioned 
above, the normalized electric field is 0.003 and 0.004, re- 
spectively. For such small normalized values, and also in 
view of the fact that the electron mean-free-path is orders 
of magnitude below the parallel wave-length, the effect of 
the acceleration on the bulk electrons is negligible. 

The electron parallel velocity in such a time-varying 



v(t) = v 



eE z 



m e (k z Vo — LOr) 



{sin[fc z 2 + (k z vo - u r )t\ 



i(k z z )} . 



(20) 



Here, vo and zo are the starting electron velocity and posi- 
tion in the parallel direction, respectively. Clearly, the ac- 
celeration of every separate particle is dependent on its par- 
ticular velocity Vq, the lowest being for those with vo — 0. A 
strong acceleration will take pla ce for resonant particles sa t- 
isfying v = uo r /k z [see also in iFletcher fc Hudson! ((2008)] . 
For the parameters used above A z = 10 km and \ y — 0.5 m, 
we have ui r — 254 Hz, and ^/uj r = 0.26. Thus, the resonant 
particles are those with Vq = 404 km/s ~ 0.1i; re . 

The acceleration of electrons with different starting ve- 
locities is seen in Fig. 7, where the achieved velocity v(t) 
(normalized to w re ) is presented for three electron popula- 
tions with vo/v T e = 0, 0.05, and 0.08. The figure describes 
the different effective oscillation periods 2n/(k z vo —uj r ) and 
the amplitudes of the three electron populations within the 
same wave period t e = 2n/u) r — 0.025 s of the electric field 
with <j> — 80 V (which coincides with the electron population 
with Vq — 0). 

The given time variation implies electron fluxes in both 
directions along the z-axis. For the line labeled HI in Fig. 7, 
we have the electron velocity Av ~ 40i>Te corresponding 
to the energy of 10~ 14 J, or around 69 keV. This oscilla- 
tory acceleration/deceleration is expected to considerably 
contribute to electron heating only if i/ e 2> uj t . The pro- 
cess is also accompanied by a radiation, with the energy 
jPanofskv fc Philips|[l978h radiated by an electron 



AS 
At" 



1 2e 2 (Av/At) 2 -(vx Av/Atf/c 2 



47T£o 3C 3 



(1 -v 2 /c 2 ) 3 



Setting as an example At — t e , for this particular case we 
obtain a very small energy AE = 1.4 ■ 10 -35 J radiated by 
an electron in one act of deceleration. 

Because the wave-field is, in fact, created by the plasma 
particles, this acceleration/deceleration should have a feed- 
back effect onto the wave itself. Clearly, it may affect the 
wave amplitude, especially for shorter wavelengths. In this 
case, the escaping electrons imply a lower amount of them 
remaining to shield the ion perturbations, and the pertur- 
bations should be increased. Yet, this all could be captured 
consistently only by numerical tools. 



4.2 Plasma heating by the drift wave 

From Eq. (|18p . it is seen that for a wave frequency below 
uj,e, the dissipation rate can formally be written as uji — 
\fei \ — \"/ion\, and two parallel mechanisms of plasma heating 
are in action here. 

The term \*yi 0n \ is responsible for the Landau dissipation 
of the wave energy and, consequently, for the heating of the 
plasma. So, as long as the density gradient is present, there 
is a continuous precipitation of energy from the wave to the 
plasma. Note that similar heating due to the Landau dissi- 



pation of the ion acousti c mode (|D'Angelolll968l ) predicts a 
stronger heating of ions <|Revathvl[l97^ r in agreement with 
observations. 
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On the other hand, the term |7 e ;| produces a growth 
of the wave, and this implies another (stochastic) heating 
mechanism that also involves single particle interaction with 
the wave. This is a process descri b ed and experimentall y 
verified in iMcChesnev et al.1 (|l987h : ISanders et ail (|l99ST ) . 
For drift wave perturbations of the form <f>{x) cos(k y y + k z z — 
ut), \k y \ \k z \, one finds the ion particle trajectory in the 
wave field from the following set of equations: 

dx/dr = T, x = k y x, T = k y y, r = Slit 

d 2 T/dr 2 = -T + [mik 2 cj>/(eBa)] sin(T - rw/ft). (21) 

The analysis from I Sanders et all (|l998T ) reveals that stochas- 
tic heating takes place for a sufficiently large wave ampli- 
tude, more precisely for 

a = k 2 y p 2 -^r>l. (22) 

The maximum achieved bulk ion velocity is shown to be 
proportional to the wave amplitude and is given by 

v max ~ [k 2 y ple(j>/{KTi) + l.9]SU/ky. (23) 

Ideally, this requires 

hel\ > \fi m \, (24) 

so that the wave amplitude may grow and at some point 
both heating mechanisms may take place, simultaneously. 
The condition (|24|l can easily be satisfied in view of the 
almost unlimited range of the parallel wave number k z , 
so that the ion Landau damping can be made small, i.e., 
\u/k z \ > v T i- 

In the stochastic heating due to the drift wave, the ions 
move in the perpendicular direction to large distances and 
feel the time- varying field of the wave due to the polarization 
drift v p = (dE/dt)/(SliB ) [the sixth term in Eq. ©], and 
as a result their motion becomes stochastic. In other words, 
the polarization drift of the ions becomes comparable to 
the E x B drift, and the displacement due to the polariza- 
tion drift is comparable to the wavelength. The single-ion 
motion becomes chaotic because of trapping (in the wave 
potential well) and de-trapping due to the magnetic field. 
The polarization drift is in the direction of the wave num- 
ber vector, which emphasizes the crucial electrostatic nature 
of the wave in the given process of heating. Also impor- 
tant to stress is that in this scenario the stochastic heating 
is highly anisotropic, and it takes place mainly in the di- 
rection normal to the magnetic field Bo (both the x— and 
j/-direction velocities are stochastic). In the same time, in 
view of the mass difference and the physical picture given 
above, this heating scenario predominantly acts on ions. The 
heating is maximal in the areas of strong density gradients 
(the areas of maximal drift wave activity), and also propor- 
tional to the strength of the magnetic field (the stronger the 
field the more localized heating). All these facts have been 
confirmed experimentally and, in the solar case, satisfy the 
observational constraints for the coronal heating mechanism 
discussed in the introduction. 

As a matter of fact, in application to coronal plasma, the 
indication or proof that the described heating really takes 
place would be: i) an ion temperature anisotropy Tt± S> T iz , 
ii) a possibly higher ion temperature in comparison to elec- 
trons, and iii) a better heating of heavier ions. Observa- 
tions show that i) may be taken rather as a rule than as 




Figure 8. The stochastically increased ion temperature T e t t = 
m i v2 nax/(3 K ) ( m millions K) in terms of the perpendicular wave- 
length X y and the ion mass. 

an exception [ s ee in iLi et al.l (|l99cf ). ICuseri et al.l (|l999h . 
ICranmer et al.l (|200Sf )]. i.e.. the perpendicular stochastic 
heating is more dominant compared to the parallel heating. 
This may also be expected from Fig. 3, where the Landau 
resonance ui/k z ~ v T i takes place for short wavelengths. For 
example, in the case of line c there, having v T i — 91 km/s, 
the ion resonance takes place at X z ~ 2.2 km, i.e., in the do- 
main where the wave is strongly damped and will not appear 
at all. There are also numerous indications that confirm the 
features ii) and i ii). A s an example, we refer to graphs from 
lHansteen et al.1 (|l997T l. where T e < T H < T He throughout 
the cor ona and the s olar w ind. Similar results may a lso be 
seen in ICuseri et al.l (|l999h and ICranmer et al.l ((2008) , and 
in references cited therein. 

A stronger heating of heavy ions can be understood 
from Eq. (|23[) and after expressing the effective temperature 
in terms of the ion mass T e ff(m,i) = rrnVm ax /(3K). From 
the derivative dT e f /(jn») /dm, > 0, we find that the heating 
increases with the ion mass if: 

fe^(^) 2 >1.9. (25) 

For X y = 0.5 m, X z = 10 km, L n — 100 m, we obtain 
from Eq. (JTHJ Ui/u T = 0.26 and u) r = 254 Hz. Note that 
in this case a ~ 1 and the stochastic heating is in action. 
Assuming small starting perturbations e<j>/(K,Ti) = 0.01, i.e., 
4> — 0.86 V, the value <f> — 60 V is achieved within r g = 
0.06 s. Those were the reasons for the amplitudes of (j> used 
in the preceding section. 

For 4> = 60 V we then have T t (X y ,p) = 0.881 + 
0.057/i/Ay + 1.78X y /fj, (normalized to the starting tempera- 
ture Ti = 10 6 K). In Fig. 8, we plot the obtained temper- 
ature in terms of the ion mass (normalized to the proton 
mass) , and the perpendicular wavelength X y . It is seen that 
for short A^, and therefore for fast growing modes, the heat- 
ing is always larger for heavier ions. 

Using the same starting set of parameter values as 
above, by using Eqs. (|22[) and (|23|l . in Table 1 we calculated 
the effective increase of the ion temperature for hydrogen 
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and helium for several wave amplitudes and perpendicular 
wavelengths. The results confirm and quantify the conclu- 
sions drawn above. The stronger heating for helium (values 
in brackets) in the short wavelength range is because the 
condition (|25[) is easily satisfied. Obviously, the proposed 
mechanism is more than efficient enough to heat the ions in 
the solar corona to the observed temperatures. 

Assuming again small perturbations e<f>/(nTi) = 0.01 
for A = 0.5 m, L n = 100 m we have the above given values 
for u)i,u) r , and the value (j> — 60 V achieved within r 9 = 
0.06 s. The maximum energy released per unit volume is 
E raal = norriiv'} nax /2 — 0.04 J/m 3 . The energy release rate 
Fmax = ^max/fg ~ 0.7 J/ (m 3 s) amounts to 4 orders of 
magnitude above the necessary value. 

However, for L n = 100 km (i.e., setting s = 1000) we 
obtain uJi = 0.07 Hz, lo t — 0.25 Hz, r g = 64 s, and conse- 
quently Tmax = 6.3 ■ 10 -4 J/(m 3 s), that is around 10 times 
the value presently accepted as necessary. Hence, the heat- 
ing rate in large magnetic loops comes close to the required 
values. Similar estimates may be made for still larger L n , yet 
the conditions under which the previous expressions are de- 
rived become violated and a numerical approach is required 
in this case. 

For the typical achieved effective temperatures from 
Table 1, we have a heating rate of ions of the order of 
10 7 — 10 s K/s, whi ch is similar to the he a ting rate obtained 
in the experiments ISanders et al.l i| 19981 ); iMcChesnev et ail 
(1987). Observe also that the magnitude of the electric field 
which we are dealing with, is of the same order as in the 
experiments. 

However, some effects that are not included here 
may reduce the energy yield, especially at shorter spatial 
scales. They r ead a s follows. In reality, the nonlinearity 
iLee fc Okudal J 19761 ) and collisions lead to the radial flat- 
tening of the density profile in the region occupied by the 
wave, resulting in the saturation of the growth. The flat- 
tening in some region Ar occupied by the wave, around a 
point ro, leads to the saturation of the instability in this 
particular region. However, the process is accompanied by a 
simultaneous steepening of the density profile outside the re- 
gion ro± Ar/2. These newly created (and also even steeper) 
density gradients will support the excitation of new modes 
now at different positions in radial direction. So here we have 
a sort of 'double cascade': the regions affected by the excita- 
tion of waves (and the heating) are shifting radially in both 
positive and negative directions. In other words, the start- 
ing instability (and heating) initiated around the position 
ro (and in the same time extraordinarily elongated in the 
axial direction because k z /k y <C 1) will have the tendency 
of spreading radially. 

The energy diffusion (due to any reason) in the perpen- 
dicular direction may reduce the local effects of the demon- 
strated strong heating rate. The earlier given Eq. (I19f) de- 
scribes the coupling with the Alfven wave, that is propor- 
tional to k z k y . All these effects (the nonlinearity, collisions, 
particle acceleration by the electric field, diffusion, and cou- 
pling with the Alfven wave) will more effectively act on 
short scales. All of them, except the particle acceleration, 
will tend to reduce the mode amplitude. The actual val- 
ues for r are thus expected to be below Y ma x. Therefore, 
the apparently too large release of energy at short scales, as 
formally obtained above, may in reality be considerably re- 



duced. Clearly, more accurate estimates and more detailed 
description may be obtained only numerically. 

Some additional features of the drift wave may also play 
a crucial role. First, the mode becomes easily nonlinear, and, 
second, within the nonlinear theory it allows a double cas- 
cade in fc-space, i.e., the transport of the wave energy both 
towards large and short wavelengths. Hence, due to non- 
linear three-wave parametric interaction, a mode growing 
for certain uo, ko will tend to excite modes at rather differ- 
ent (both smaller and larger) ujj , kj , and in the end those 
will always include the modes that heat the system most 
effectively. For example, the shape of Eq. (I23p reveals that 
v m ax ~ a\k v + (12/ky and, in terms of k y , v ma x has a min- 
imum at k'i — 1.9KTi/(pie(j>)- Hence, even if a mode with 
wave number k c does not produce a so large maximum bulk 
velocity v ma x, it will nonlinearly excite modes with k y far 
from k c and there will be a stronger heating with this non- 
linearly generated mode. 



4.3 Heating of cool corona 

In sections 4.1 and 4.2, we have demonstrated the possibility 
for sustainable coronal temperatures of around million K 
with the discussed drift wave mechanism. Assume now a 
relatively cool starting corona, with a temperature of only 
10 4 K (the value taken high enough in order to neglect the 
presence of neutrals) and let us calculate the time that is 
needed to achieve the temperature of 1 million K, using the 
same mechanism. 

We first chose the least favorable value k y — k c calcu- 
lated above. This yields a = 1.9 for any mode amplitude <f>, 
so we have the stochastic heating in action. For these a, k c , 
from Eq. ((23j) we obtain v max = S.8D. l p l [e ( f>/(1.9KT i )] 1/2 . 
The effective temperature of 1 million K implies that the 
mode amplitude has the value (j> = 4> m = 34 V. This then 
yields k y = k c = 23.14 1/m and A c = 0.27 m. 

Hence, choosing A = A c , and A z = 50 km, L n = 100 m, 
and n = 10 15 m" 3 , from Eqs. (JTTJ) and {THJ we find uo r — 
18.3 Hz and 7 = 0.75 Hz. Assuming a small starting pertur- 
bation of the potential <j> = 0.0086 V, i.e., e<f>/(K,Ti) = 0.01, 
with the given growth rate we find out that the mode am- 
plitude <j) m (at which the temperature achieves the required 
value of million K due to the stochastic heating), is reached 
at longest within lis. For any other value of \ y we will have 
a shorter growing and heating time. 



5 KINETIC CURRENT-DRIVEN DRIFT 
WAVE INSTABILITY 

In the presence of an additional energy source, in the form of 
an electron current in the direction of the magnetic field vec- 
tor, the drift wave becomes even more unstable as compared 
to the previous collisionless kinetic instability. The mode fre- 
quency and the growth rate, with so me additional features 
and details (Hatakc vama et al.lll980l ). are then given by: 



A (&i 



[1 



1 + [1-A (&;)]»- 
1 + [1 - A (&i)]?7_L + Ao(&i)77„ kin 



2 

Tiz 



l + [l-Ao(6i)]»7i 



2uji 



(26) 



The universally growing mode in the solar atmosphere: coronal heating by drift waves 11 



Table 1. Plasma heating for hydrogen and helium (in 
brackets) ions for several values of the perpendicular 
wavelength and the wave amplitude <j>. The maximum 
stochastic velocity is given by Eq. (1 23 I t and the effec- 
tive temperature obtained by heating is m"»jL ax /(3re). 



<t> = 60 [V] 


\ y [m] 


Vmax [m/s] 








T*ff [ K ] 


0.3 


2.12 • 10 5 (1.50 


10 5 ) 


1.82 


10 6 


(3.50 • 10 6 ) 


0.5 


2.20 ■ 10 5 (1.10 


10 5 ) 


1.96 


10 6 (2.01 • 10 6 ) 


0.8 


2.79 ■ 10 5 (1.78 


10 5 ) 
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10 6 


(1.78 • 10 6 ) 
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Vmax [m/s] 
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Teff [K] 
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(5.78 • 10 7 ) 
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2.45 ■ 10 5 (1.37 
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2.43 
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(3.02 • 10 6 ) 
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10 6 
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Teff [K] 
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10 6 
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2.71 ■ 10 5 (1.62 
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10 6 


(4.23 • 10 6 ) 
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3.10 • 10 5 (1.36 
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rV2 



1 + [1- A (&i)]77j 



1 + 



7] z A k^VTiz/^*ep 



— + — I 

. VTe K z VTe 



-Ao(&i 



2[l + [l-Ao(6i)]»7i 



r 2 

x exp|— a>* e 



,/(k 2 z VTiz)]} 



(27) 



Here, uo is the electron current, the starting ion tem- 
peratures in the perpendicular and parallel direction are 
allowed to be different , and we use the notation from 
lHatakevama et al] (|l980) n±, z = T e /Ti,±, x , v% iz — uTiz/rrii, 
and oj* ep = cj* e A (fe;)/[l + (1 _ Ao(6i))r/x]- 

In solar corona environment, such electron currents may 
appear in the processes of magnetic reconnection. An ex- 
ample with the s ame phenomenon for th e terrestrial atmo- 
sphere is given in lSerizawa fc Satol (|l984T ). 

We calculate the growth rate (|27[) numerically by tak- 
ing B = 10" 3 T, n = 10 13 m" 3 , L n = [{dn / dx) /no]- 1 = 
1 km, T e = T ix = 10 6 = 2T iz , and A 2 = 200 km. This differ- 
ent set of parameters is chosen only for a change, in view of 
the comments given earlier in Sec. 4 (see also Fig. 6). The re- 
sult is presented in Fig. 9, where the growth rate u>i(X y ,uo) 
is normalized to the wave frequency oj r (A H ,Mo). We observe 
that for the given set of parameters, the current addition- 
ally increases the growth rate. This is seen by comparing the 
limits uo = and uq = 10, for e.g. X y = 15 m, where the 
growth rate is increased by about a factor 20. However, the 
effect is of less importance for shorter wavelength (see the 
limit X y = 5 m). In the given example, the mode frequency 
is below w* e , so that the wave is growing anyhow, and this is 
due to the kinetic effects presented in the previous sections. 

The real importance of such an electron flow is com- 
pletely different in the limit to ~ oj* e when the instability 
sets in provided that the electron current exceeds a certain 
threshold. This may be checked for example by setting larger 
values for Bo when the kinetic instability from the previous 




Figure 9. Growth rate LUi(X y , «o)/w(A s , no) from Eqs. J26D and 
J27} for the kinetic + current-driven drift-wave instability in the 
case B = 0.001 T. 



sections vanishes (because for the given parameters w* e is re- 
duced so that lo ^ w* e ). Yet, the wave instability reappears 
for a sufficiently large wo, but it is now the current-driven 
one. Ex actly such a sort of b ehavi or was verified experimen- 
tally in lHatakevama et all (| 1980T) . Hence, for example we 
take the magnetic field for one order of magnitude larger 
Bo = 10~ 2 T, and the result is presented in Fig. 10. The 
contour lines of oJi(X y , uo)/u(X y , uo) in Fig. 10 yield only 
the current-driven drift wave for sufficiently large X y . 
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£ 1.1 



X [m] 

y 




Figure 10. Contour lines of ^(A^, uq)/w(\ v , uq) from Eqs. J26D 
and l|27|l for Bo = 0.01 T, showing the current-driven drift-wave 
instability in the vicinity of threshold. 



5.1 Plasma heating 

The heating of the plasma by this curre nt-driven drift-wave 
instab ility, as experimentally measured in lHatakevama et alj 
(1980), has similar properties as the one discussed in the pre- 
vious text. This implies the following facts: i) the anisotropy 
in the ion heating (Tjj. > Ti z , ii) the heating time being 
comparable to the growth time of the drift wave and pro- 
portional to the wave energy (square amplitude), iii) the 
perpendicular heating being associated with the leading or- 
der perturbed perpendicular drift, while the parallel heating 
is associated with the ion Landau damping. 

We stress that the colhsional counterpart of the above 
described current-driven drift-wave instabilit y, applicable to 
the lo wer solar atmosphere, may be found in Ellis Motley! 
( 1978) where it has been explained and also experimentally 
verified. 



6 ION-CYCLOTRON DRIFT WAVE 

In the case of high-frequency electrostatic perturbations 
satisfying the condition uo r ~ Qi, and in the presence of 
a density gradient, we have an ion-cyclotron drift (ICD) 
wave, also calle d the ion Bernstein mode in the literature 
Ichim aru 1980) . This mode has b een theoretically predicted 



Figure 11. Contour plot of S(k y ,L n ) from Eq. J28). The ICD 

mode is growing for S(k y ,L n ) > 0. 



■Mikhailovskii fc Timofeevl 1 19631 ) and experimentally veri- 
fied (|Hendel fc Yamaddll97lT lo ng ago. The kinetic descrip- 



tion given in llchimarul (Il98d); [Mikhailovskii fc Timofeevl 



|l963MHendel fc Yamadal (|l974h reveals the same instabil- 
ity mechanism as in the case of the kinetic drift wave insta- 
bility presented before. However, in order for the instability 
to take place, the drift-wave branch and the ion-cyclotron 
branch must get close to each other. In that case, instead 
of the intersection of two dispersion lines, we have only one 
complex-conjugate solution. Clearly, to have this, the equi- 
librium density scale length must be very short to make the 
frequency of the drift part high, this because for the drift 
wave u> r ~ 1/L n . An application to solar plasma of this in- 



st ability is d i scuss ed in lVranies fc Poedtsl |2008a). As shown 
m llchimarul (|l980h . the ICD mode grows if the following in- 
stability condition is satisfied: 

S(k y ,L n ) = 2ky\dU!l(LJl +U)2)5 

~ {uj 1 ~uj. 2 ) 2 {l + k 2 v \ 2 d -28) >0. (28) 

Here, Wi = w„ e /(l+fc^), lo 2 = Oi [l + 5/(1 + k 2 y \ 2 d )] , = 
p 2 m e /m i + \ 2 de Ti/T e , \ de = v Te /uj pe , S = l/[(2n) 1/2 k y pi] 

We stress that the ion-cyclotron mode itself has been 
discussed a lot in the past in the context of p roblems related 
to the heating of the so l ar corona, see e.g. ICranmer et all 
l|l999l) . ICranmerl l|2000l ). iHollweg fc Isenberd |2002l ). The 
reasons for this are the same as those discussed in the pre- 
vious sections: the evidence obtained from in situ measure- 
ments in the solar wind and coronal holes of resonant ion 
cyclotron heating, and a preferential heating of coronal ions 
(with respect to electrons), that in the same time is most 
dominant in the direction perpendicular to the magnetic 
field lines. In this context, the damping of such IC waves 
is believed to be a good candidate for th e consequent coro - 
nal heating and solar wind acceleration l|Markovskiill200ll ). 
Yet, as in many wave-heating scenarios of the corona pro- 
posed in the past, there is the problem of the source for 
the required generation of such IC waves throughout the 
corona. The effects proposed a s sources for the IC mode 
are global resonant MHD mo des fMarkovskii 2001), currents 
l|Forslun J200ll : |Toichilll97ll ) , etc. However, as a rule, these 
effects themselves need some source, so the problem is not 
solved but merely shifted to ano ther problem. 

The analysis performed in IVranies fc Poedts! (|2008al ) 
and in the previous sections, demonstrates that such a source 
(the density gradient) exists, and the mechanism which it 
implies is well known in the literature but not used or stud- 
ied in the context of the solar plasma, and it is able to pro- 
duce growing ICD modes on massive scales. For example, in 
Fig. 9 we give the contour plot of S(k y ,L n ) from Eq. (f28l) . 
The figure shows a narrow range of L n and X y within which 
the ICD is growing, S(k y ,L n ) > 0. The parameters arc 
the same as before: Bo = 10 -3 T, no = 10 13 m -3 , and 
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T e — Ti — 10 6 K. The instability shown here implies short- 
scale density inhomogeneities. However, this scale can be 
increased for example by decreasing the magnetic field in- 
tensity. Setting 10 -4 T and also reducing the number den- 
sity by one order of magnitude, yields the necessary density 
scale length L n for the unstable modes of the order of 60 
meters. Observe that such reduced values may be obtained 
e.g. by moving further along a flux tube, implying that the 
mechanism may work at various altitudes in the solar atmo- 
sphere. 

6.1 Stochastic heating of plasma by ICD wave 

A stochastic heating mechanism, similar to the case dis- 
cussed previously for the drift wave, is also known to exist in 
th is particular frequency domain. This has been discussed 
in iKarnev fc Bers I l| 19771 ). It is shown that the stochastic 
heating by a perpendicularly propagating electrostatic wave, 
satisfying the condition u) T ~ Ifli, I = 1, 2, • • •, takes place if 
the mode amplitude exceeds the limit: 

E>?fp) 1/3 . (29) 

4 kj_ \LO r J 

The condition is obtained after analyzing the same equations 
as Eq. (f2"Tj) . This yields t he wave potential <j> = j? fM 2/3 /4 = 
24/ 2/3 . The results from IVranies fc Poedtsl ()2008al ) provide 
the maximum mode growth rate uJi/uj r — 0.1. For the 
given magnetic field we have fli ~ 10 Hz for protons. 
Now, assuming small starting electrostatic perturbations 
ecf)/(hiTi) ~ 0.01, and setting I = 1, we find out that the 
required wave potential amplitude is achieved within about 
3 ■ 10 -4 s only! Hence, in the given coronal environment 
the stochastic heating condition (|29|) is practically instan- 
taneous. It develops in the way already described in the 
pervious text. 



7 CONCLUDING COMMENTS 

We have presented a few (out of many) physical effects that 
make the drift mode growing. There exist many more phe- 
nomena that lead to the growth of drift waves but these 
are not discussed here at all, e.g. the temperature gradi- 
ent driven drift wave instability, shear flow instability, etc. 
As a matter of fact, the presence of a temperature gra- 
dient implies an additional source of energy and an in- 
stability provided that th e temperature and density gra- 
dients have opposite sign l|Rudakov fc Sagdeev I 1961). i.e., 
d\nT/dlnno < 0. One particularly strong mechanism lead- 
ing to an instability is a sheared flow. The term 'shear flow' 
refers to the flow of the plasma as a whole along the mag- 
netic field vector and having in the same time a gradient 
in the perpendicular direction. Such a current-less insta- 
bility [known also as D'Angelo mode (|D'Angelolll965l )] has 
an extra energy source in the flow gradient. A kinetic the- 
ory a nalysis of this instability as given in IGanguli et al. I 
l|l994h . whil e its most rece n t exp erim ental verification may 
be found in iKaneko et al. I (I2003T ). In ISaleem et al. I (|2007l ) 
such a shear flow instability is discussed within the fluid 
theory in application to coronal spicules, showing strongly 
growing modes with the possibility for uji > u r . Because 
of the radial density gradients, solar spicules are an ideal 



nursery for the drift wave instability. Growing drift waves 
are simply carried upwards by the plasma flow inside each 
of them. In the same time they are numerous (at least 10 5 
of them are present throughout the Sun at every moment). 
There are also evidences of the presence of plasm a flows 
along the magnetic flux loops ISchriiver et al. | (|1999| ). While 
the magnitude of these axial flows is determined (showing 
subsonic flows with velocities up to 100 km/s), their even- 
tual inhomogeneity in the radial direction remains an open 
question. Yet, there are no obvious physical reasons that 
would exclude them. The inflow of more dense plasma along 
a magnetic flux tube from lower (cooler) layers may then im- 
ply the simultaneous presence of a density gradient (towards 
the axis) and a temperature gradient (outward oriented) . In 
such a geometry, the mentioned temperature-gradient in- 
stability may develop, or/and it may be combined with the 
shear-flow instability. 

In reality, and in particular in the heating processes in 
the solar corona, the simple growth-damping process of a co- 
herent wave, and the consequent heating due to the resonant 
interaction with plasma particles, can not be a completely 
accurate picture. It is rather an interplay between various 
processes that may happen in the same time and place, like 
collisions, mode growth, kinetic effects, nonlinearity and tur- 
bulence, and stochastic ion and electron heating. The results 
presented here show that the stochastic heating related to 
the drift wave is a powerful mechanism that switches on for 
sufficiently large mode amplitudes. The properties of the re- 
sulting heating process are such that they are consistent with 
all the observed features of the plasma heating in the solar 
atmosphere. On the other hand, there is plenty of energy for 
the mode itself, stored in the omnipresent density gradients. 
To our knowledge, there is no other available heating model 
like the one presented in this work, that is so clearly able to 
fulfill all these requirements. 

The presented mechanism removes the necessity for 
explaining the most crucial problems of coronal heating: 
namely i) how the (right amount of) energy is transmit- 
ted from 'sources below' the photosphere to the corona, and 
ii) how this energy is dissipated locally in the corona (at 
the right rate). Instead, we showed, first, that a sufficient 
amount of energy for driving drift modes is already present 
in the corona, and, second, that it is naturally transmitted 
to the different plasma species by well known effects that 
are, however, beyond the standardly used models and the- 
ories. Hence, the proposed mechanism is based on a novel 
paradigm that allows a self-consistent solution model. The 
heating mechanism implies instabilities on time and spatial 
scales that are currently not directly observable by space 
probes. However, all the effects presented here are directly 
experimentally verified under laboratory conditions. Their 
indirect confirmation in the context of the solar corona seems 
to be also beyond doubts. This is because the consequences 
of the heating process, as enlisted earlier in the text (temper- 
ature anisotropy, better heating of heavier ions, hotter ions 
than electrons), are indeed verified by satellite observations. 



ACKNOWLEDGMENTS 

The results presented here are obtained in the framework 
of the projects G.0304.07 (FWO-Vlaanderen) , C 90347 



14 J. Vranjes and S. Poedts 



(Prodex), GOA/2009-009 (K.U.Leuven). Financial support 
by the European Commission through the SOLAIRE Net- 
work (MTRN-CT-2006-035484) is gratefully acknowledged. 
TRACE (Fig. 1) is a mission of the Stanford-Lockheed In- 
stitute for Space Research, and part of the NASA Small 
Explorer program. 



REFERENCES 

Aschwanden, M. J., 2001, ApJ, 560, 1035 

Aschwanden, M. J., 2004, Physics of the Solar Corona. 

Springer- Verlag, Belin-Heidelberg-New York 
Aschwanden, M. J., Wiilser, J. P., Nitta, N. V., Lemen, R. 

R., 2008a, ApJ, 679, 827 
Aschwanden, M. J., Lee, J. K., Gary, G. A., Smith, M., 

Inhester, B., 2008b, Sol. Phys., 248, 359 
Bittencourt, J. A., 1995, Fundamentals of Plasma Physics. 

Sao Jose dos Campos, Brazil 
Brossier, P., Deschamps, P., Gravier, R., Pellat, R., Re- 

naud, C, 1971, Phys. Rev. Lett., 26, 124 
Chen, F. F., 1988, Introduction to Plasma Physics and 

Controlled Fusion. Plenum Press, New York 
Cranmer, S. R., Field, G. B., Kohl, J. L., 1999, ApJ, 518, 

937 

Cranmer, S. R., 2000, ApJ, 532, 1197 

Cranmer, S. R., Panasyuk, A. V., Kohl, J. L., 2008, ApJ, 
678, 1480 

Cuseri, I., Mullan, D., Poletto, G., 1999, Space Sci. Rev., 
87, 153 

D'Angelo, N., 1965, Phys. Fluids, 8, 1748 
D'Angelo, N., 1968, ApJ, 154, 401 

De Pontieu, B., Mcintosh, S. W., Carlsson, M., Hansteen, 
V. H., Tarbell, T. D, Schrijver, C. J., Title, A. M., Shine, 
R. A., Tsuneta, S., Katsukawa, Y., Ichimoto, K., Sue- 
matsu, Y., Shimizu, T., Nagata, S., 2007, Science, 318, 
1574 

Dreicer, H., 1959, Phys. Rev., 115, 238 

Ellis, R. F., Motley, R. W., 1978, Phys. Fluids, 17, 582 

Fletcher, L., Hudson, H. S., 2008, ApJ, 675, 1645 

Forslund, D. W., 1970, J. Geophys. Res., 75, 17 

Ganguli, G., Keskinen, M. J., Romero, H., Heelis, R., 

Moore, T., Pollock, C, 1994, J. Geophys. Res., 99, 8873 
Hansteen, V. N., Leer, E., Holtzer, T. E., 1997, ApJ, 482, 

498 

Hasegawa, A., 1975, Plasma Instabilities and Nonlinear Ef- 
fects. Springer- Verlag, Berlin Heidelberg 

Hatakeyama, R., Oertl, M., Mark, E., Schrittwieser, R., 
1980, Phys. Fluids, 23, 1774 

Hendel, H. W., Yamada, M., 1974, Phys. Rev. Lett., 33, 
1076 

Hollweg, J. V., Isenberg, P. A., 2002, J. Geophys. Res., 107, 
10.1029 

Ichimaru, S., 1980, Basic Principles of Plasma Physics. 
The Benjamin/Cummings Publish. Comp., Reading, Mas- 
sachusetts 

Kadomtsev, B. B., 1965, Plasma Turbulence. Acaddemic 
Press, London 

Kaneko, T., Tsunoyama, T., Hatakeyama, R., 2003, Phys. 

Rev. Lett., 90, 125001 
Karney, C. F. F., Bers, A., 1977, Phys. Rev. Lett., 39, 550 
Karovska, M., Habbal, S. R., 1991, ApJ, 371, 402 



Klimtchuk, J. A., 2006, Solar Phys., 234, 41 
Lee, W. W., Okuda, H., 1976, Phys. Rev. Lett., 36, 870 
Li, X., Habbal, S. R., Kohl, J. L., Noci, G., 1998, ApJ, 501, 
L133 

Markovskii, S. A., 2001, ApJ, 557, 337 

McChesney, J. M., Stern, R. A., Bellan, P. M., 1987, Phys. 

Rev. Lett., 59, 1436 
Mecheri, R., Marsch, E., 2008, A&A, 481, 853 
Mikhailovskii, A. B., Timofeev, A. V., 1963, Sov. Phys. 

JETP, 17, 626 

Milligan, R. O., Gallagher, P. T., Mathioudakis, M., 
Keenan, F. P., Bloomfield, D. S., 2005, MNRAS, 363, 259 

Narain, U., Ulmschneider, P., 1990, Space Sci. Rev., 54, 
377 

Nishida, Y., Ishii, K., 1974, Phys. Rev. Lett., 33, 352 
November, L. J., Koutchmy, S., 1977, ApJ, 466, 512 
Pandey, B. P., Wardle, M., 2008, MNRAS, 385, 2269 
Panofsky, W. K. H., Philips, M., 1978, Classical Electricity 

and Magnetism. Addison- Wesley, Reading 
Regnier, S., Parnell, C. E., Haynes, A. L., 2008, A&A, 484, 

L47 

Revathy, P., 1977, Sol. Phys., 53, 445 

Rudakov, L. I., Sagdeev, R. Z., 1961, Sov. Phys. Doklady, 
6, 415 

Sanders, S. J., Bellan, P. M., Stern, R. A., 1998, Phys. 

Plasmas, 5, 716 
Saleem, H., Vranjes, J., Poedts, S., 2007, A&A, 471, 289 
Schrijver, C. J., Title, A. M., Berger, T. E., Fletcher, L., 
Hurlburt, N. E., Nightingale, R. W., Shine, R. A., Tarbell, 
T. D., Wolfson, J., Golub, L., Bookbinder, J. A., Deluca, 
E. E., Mcmullen, R. A., Warren, H. P., Kankelborg, C. 
C, Handy, B. N., De Pontieu, B., 1999, Sol. Phys., 187, 
261 

Serizawa, Y., Sato, T., 1984, Geophys. Res. Lett., 11, 595 

Suzuki, T. K., 2004, MNRAS, 349, 1227 

Toichi, T., 1971, Sol. Phys., 18, 150 

Vranjes, J., Poedts, S., 2004, Phys. Plasmas, 11, 891 

Vranjes, J., Poedts, S., 2006, A&A, 458, 635 

Vranjes, J., Poedts, S., 2006, Phys. Plasmas, 13, 122103 

Vranjes, J., Poedts, S., 2007, Phys. Plasmas 14, 112106 

Vranjes, J., Pandey, B. P., Poedts, S., 2007, Phys. Rev. 

Lett., 98, 049501 
Vranjes, J., Poedts, S., 2008a, A&A, 482, 653 
Vranjes, J., Poedts, S., 2008b, Phys. Plasmas, 15, 034504 
Vranjes, J., Poedts, S., Pandey, B. P., De Pontieu, B., 

2008a, A&A, 478, 553 
Vranjes, J., Kono, M., Poedts, S., Tanaka, M. Y., 2008b, 

Phys. Plasmas, 15, 092107 
Vranjes, J., Poedts, S., 2009, Phys. Plasmas, 16, 022101 
Warren, H. P., Ugarte-Urra, I., Doschek, G. A., Brooks, D. 

H., Williams, D. R., 2008, ApJ, 686, L131 
Weiland, J., 2000, Collective Modes in Inhomogeneous 

Plasmas. Institute of Physics Pub., Bristol 
Woo, R., Habbal, S. R., 1997, ApJ, 474, L139 
Woo, R., 1996, Nature, 379, 321 

This paper has been typeset from a TgX/ DTj^X file prepared 
by the author. 



1-Or 




A z /s [km] 



